Direct measurement of the Wigner function by photon counting 
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We report a direct measurement of the Wigner func- 
tion characterizing the quantum state of a light mode. The 
experimental scheme is based on the representation of the 
Wigner function as an expectation value of a displaced pho- 
ton number parity operator. This allowed us to scan the phase 
space point-by-point, and obtain the complete Wigner func- 
tion without using any numerical reconstruction algorithms. 

PACS Number(s): 42.50.Ar, 03.65.Bz 



Among many representations of the quantum state, 
the Wigner function offers an appealing possibility to de- 
scribe quantum phenomena using the classical-like con- 
cept of phase space 0. The Wigner function provides 
complete information on the state of a system, and it al- 
lows one to evaluate any quantum observable by phase 
space integration with an appropriate Wigner- Weyl or- 
dered expression. 

Recently, the Wigner function has gained experimen- 
tal significance due to the development of the optical 
homodyne tomography, a beautiful technique for mea- 
suring the quantum state of light pioneered by Smithcy 
et al. H and further applied by Breitenbach et al. || 
In this method, rooted in the domain of image process- 
ing, the Wigner function is a natural representation of 
the quantum state reconstructed from experimental data. 
However, the route from raw experimental results to the 
Wigner function is not straightforward. First, a sample 
of homodyne events is collected and stored. Statistics of 
these events for a fixed local oscillator phase is described 
by a marginal projection of the Wigner function. In or- 
der to retrieve the complete Wigner function, a family of 
homodyne statistics measured for a sufficiently dense set 
of local oscillator phases has to be processed using the 
sophisticated filtered back-projection algorithm. 

In this Communication we report a direct measure- 
ment of the Wigner function of a light mode. This tech- 
nique, based on photon counting, avoids the detour via 
complex numerical reconstruction algorithms. The prin- 
ciple of our measurement is entirely different from op- 
tical homodyne tomography. The Wigner function at a 
given phase space point is itself a well defined quantum 
observable [Q. Furthermore, the measurement of this 
observable can be implemented for optical fields using 
an arrangement employing an auxiliary coherent probe 
beam The amplitude and the phase of the probe 



field define the point in the phase space at which the 
Wigner function is measured. This allowed us to scan 
the phase space point-by-point, simply by changing the 
parameters of the probe field. A variation of this idea 
has been applied by Leibfried et al. JIJ to determine the 
vibrational state of a trapped ion. Here we present an ex- 
periment, which to the best of our knowledge is the first 
direct measurement of the Wigner function for optical 
fields. 

Our experiment is based on the representation of the 
Wigner function at a complex phase space point denoted 
by a as the expectation value of the following operator: 
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where D(a) is the displacement operator and \n) denote 
Fock states, n\n) = n\n). Thus, W(a) has two eigen- 
values: 2/tt and —2/tt, corresponding to degenerate sub- 
spaces spanned respectively by even and odd displaced 
Fock states. Practical means to translate this formula 
into an optical arrangement are quite simple [^[(6). The 
displacement transformation can be realized by superpos- 
ing the measured field at a low-reflection beam splitter 
with a strong coherent probe beam. The value of the 
displacement a is equal in this setup to the reflected am- 
plitude of the probe field. Furthermore, the projections 
on Fock states can be obtained by photon counting as- 
suming unit quantum efficiency. These two procedures, 
combined together, provide a practical way to measure 
the Wigner function at an arbitrarily selected phase space 
point a. 

The experimental setup we used to measure the 
Wigner function is shown schematically in Fig. |l|. In prin- 
ciple, it is a Mach-Zender interferometric scheme with the 
beams in two arms of the interferometer serving as the 
signal and the probe fields. An attenuated, linearly po- 
larized (in the plane of Fig. ^) 632.8 nm beam from a 
frequency-stabilized single-mode He:Ne laser is divided 
by a low-reflection beam splitter BS1. The weak re- 
flected beam is used to generate the signal field whose 
Wigner function will be measured. The state preparation 
stage consists of a neutral density filter ND and a mir- 
ror mounted on a piezoelectric translator PZT. With this 
arrangement, we are able to create pure coherent states 
with variable phase as well as their incoherent mixtures. 
Though these states do not exhibit nonclassical proper- 
ties, they constitute a nontrivial family to demonstrate 
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the principle of the method, which provides a complete 
characterization of both quantum and classical field fluc- 
tuations. 

The strong beam leaving the beam splitter BS1 plays 
the role of the probe field with which we perform the 
displacement transformation D(a). In order to scan the 
phase space one should be able to set freely its ampli- 
tude and phase, which define respectively the radial and 
the angular coordinates in the phase space. The am- 
plitude modulation is achieved with a half-wave plate, 
a longitudinal Pockels cell EOM1, and a polarizer ori- 
ented parallel to the initial direction of polarization. The 
phase modulation is done with the help of an ADP crystal 
electrooptic phase modulator EOM2 on the signal field. 
This is completely equivalent to modulating the probe 
field phase, but more convenient for technical reasons: in 
this arrangement optical paths in both the arms of the 
Mach-Zender interferometer are approximately the same, 
and better overlap of the signal and the probe modes is 
achieved at the output of the interferometer. 

The signal and the probe fields are interfered at a 
nearly completely transmitting beam splitter BS2 with 
the power transmission T = 98.6%. In this regime, the 
transmitted signal field effectively undergoes the required 
displacement transformation. Spurious reflections that 
accompany the beam leaving the interferometer are re- 
moved using the aperture A. Finally, the transmitted 
signal is focused on an EG&G photon counting mod- 
ule SPCM-AQ-CD2749, whose photosensitive element is 
a silicon avalanche diode operated in the Geiger regime. 
The overall quantum efficiency of the module specified by 
the manufacturer is rj > 70%. The count rate is kept low 
in the experiment, and thus the chance of two or more 
photons triggering a single avalanche signal is very small, 
and the probability of another photon arriving during 
the detector dead time can be neglected. Under these 
assumptions, each pulse generated by the module cor- 
responds to the detection of a single photon ||. The 
pulses are acquired by a computer, which also controls 
the voltages applied to the electrooptic modulators. The 
interference visibility in our setup has been measured to 
be v > 98.5%, and the phase difference between the two 
arms was stable up to few percent over times of the order 
of ten minutes. 

In Fig. H we depict the measured Wigner functions of 
the vacuum, a weak coherent state, and a phase diffused 
coherent state. The phase fluctuations were obtained 
by applying a 400 Hz sine waveform to the piezoelectric 
translator. For all the plots, the phase space was scanned 
on a grid defined by 20 amplitudes and 40 phases. The 
scaling of the radial coordinate is obtained from the av- 
erage number of photons n vac detected for the blocked 

signal path. Thus the graphs are parameterized with 
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the complex variable (3 = e l¥, n V ac , where (p is the phase 
shift generated by the phase modulator EOM2. At each 
selected point of the phase space, the photocount statis- 
tics p n {(3) was determined from a sequence of N = 8000 



counting intervals, each r = 30/xs long. The duration 
of the counting interval t defines the temporal envelope 
of the measured mode. The count statistics was used to 
evaluate the alternating sum 
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which, up to the normalization factor 2/ir is equal to 
the Wigner function of the measured state. Statistical 
variance of this result can be estimated by Var[n(/3)] = 
{1 - [U{P)] 2 }/N §. Thus, the statistical error of 
our measurement reaches its maximum value, equal to 
l/N 1 / 2 w 1.1%, when the value of the Wi gner function 
is close to zero. 

The Wigner functions of the vacuum and of the co- 
herent state are Gaussians centered at the average com- 
plex amplitude of the field, and their widths characterize 
quantum fluctuations. It can be noticed that the mea- 
sured Wigner function of the coherent state is slightly 
lower than that of the vacuum state. In the following, 
when discussing experimental imperfections, we shall ex- 
plain this as a result of non-unit interference visibility. In 
the plot of the Wigner function of the phase diffused co- 
herent state, one can clearly distinguish two outer peaks 
corresponding to the turning points of the harmonically 
modulated phase. 

There are several experimental factors whose impact 
on the result of the measurement needs to be analyzed. 
First, there are losses of the signal field resulting from 
two main sources: the reflection from the beam splitter 
BS2 and, what is more important, imperfect photodetec- 
tion characterized by the quantum efficiency 77. Analysis 
of these losses shows [|[||, that in such a case the al- 
ternating series evaluated from photocount statistics is 
proportional to a generalized, s-ordered quasidistribu- 
tion function VF(o;; s), with the ordering parameter equal 
s = — (1 — r)T)/r]T. 

In addition, the two modes interfered at the beam split- 
ter BS2 are never matched perfectly. The effects of the 
mode mismatch can be discussed most thoroughly within 
the multimode theory ]l0| ]. Here, due to limited space, 
we shall present the main conclusions and briefly sketch 
the reasoning. Let us consider the normalized mode func- 
tions describing the transmitted signal field and the re- 
flected probe field. The squared overlap £ of these two 
mode functions can be related to the interference visibil- 
ity v as £ = vj (2 — v). In order to describe the effects of 
the mode mismatch, we will decompose the probe mode 
function into a part that precisely overlaps with the sig- 
nal, and the orthogonal remainder. The amplitude of the 
probe field effectively interfering with the signal is thus 
multiplied by ^ 1 ^ 2 , and the remaining part of the probe 
field contributes to independent Poissonian counts with 
the average number of detected photons equal (1 — £)|/3| 2 . 
Consequently, the full count statistics is given by a convo- 
lution of the statistics generated by the interfering fields, 
and the Poissonian statistics of mismatched photons. A 



2 



simple calculation shows, that the alternating sum eval- 
uated from such a convolution can be represented as a 
product of the contributions corresponding to the two 
components of the probe field: 



n(/3) = cxp[-2(l-OI/?| 2 ] 



(3) 



Thus, in a realistic case II(/3) represents a Gaussian cen- 
tered at the attenuated amplitude \f£jfTctQ, and the 
width remains unchanged. This Gaussian function in 
multiplied by the constant factor exp[— 2(1 — £)?7T|ao| 2 ]- 
For our measurement, £ w 97% and r)T\ao\ 2 « 1.34, 
which gives the value of this factor equal 0.92. This result 
agrees with the height of the experimentally measured 
Wigncr function of a coherent state. 

We shall conclude this Communication with a brief 
comparison of the demonstrated direct method for mea- 
suring the Wigner function with the optical homodyne 
tomography approach. An important parameter in the 
experimental quantum state reconstruction is the detec- 
tion efficiency. Currently, higher values of this parameter 
can be achieved in the homodyne technique, which de- 
tects quantum fluctuations as a difference between two 
rather intense fields. Such fields can be efficiently con- 
verted into photocurrent signals with the help of p-i-n 
diodes. It should be also noted that an avalanche pho- 
todiode is not capable of resolving the number of simul- 
taneously absorbed photons, and that it delivers a signal 
proportional to the light intensity only in the regime de- 
scribed in this paper. However, continuous progress in 
single photon detection technology gives hope to over- 
come current limitations of photon counting fll|| . Alter- 
natively, the displacement transformation implemented 
in the photon counting technique can be combined with 
efficient random phase homodyne detection. This yields 
the recently proposed scheme for cascaded homodyning 

The simplicity of the relation (|1J) linking the count 
statistics with quasidistribution functions allows one to 
determine the Wigner function at a given point from a 
relatively small sample of experimental data. This fea- 
ture becomes particularly advantageous, when we con- 
sider detection of multimode light. Optical homodyne to- 
mography requires substantial numerical effort to recon- 
struct the multimode Wigner function. In contrast, the 
photon counting method has a very elegant generaliza- 
tion to the multimode case: after applying the displace- 
ment to each of the involved modes, the Wigner function 



at the selected point is simply given by the average par- 
ity of the total number of detected photons. Moreover, 
the dichotomic outcome of such a measurement provides 
a novel way of testing quantum nonlocality exhibited by 
correlated states of optical radiation jb| . 
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Here on the right-hand side we have made use of the 
theoretical results for imperfect detection obtained in 
Refs. [||J6). Specializing this result to a coherent state 
\(Xq) with the amplitude oto, yields: 

II(/3) =exp[-2|/3- v ^Ta | 2 -2(l-0^|ao| 2 ]. (4) [1] 
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FIG. 1. The experimental setup for measuring the Wigner function. BS1 and BS2 are quartz plates serving as low- reflect ion 
beam splitters. The quantum state is prepared using the neutral density filter ND and a mirror mounted an a piezoelectric 
translator PZT. The electrooptic modulators EOMI and EOM2 control respectively the amplitude and the phase of the point 
at which the Wigner function is measured. The signal field, after removing spurious reflections using the aperture A, is focused 
on a single photon counting module SPCM. 
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